Microscopic quark study of the 77 and rf masses 
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We show that it is necessary to go beyond the BCS (rainbow-ladder) approximation to split the 
rj and rj' masses from the tt and K masses. We determine the self-consistent set of one-quark- 
loop diagrams both for the Schwinger-Dyson quark mass gap equation and for the Bethe-Salpeter 
quark-antiquark boundstate equation. We identify the dominant diagrams, and we focus on the 
boundstate equation. We detail the Bethe-Salpeter equation, adding the dominant new diagram to 
the BCS kernel. The relevant numerical techniques are also discussed. The ideal cases of one, two 
and three light flavors, relevant to lattice QCD are also explored, together with the case of realistic 
current quark masses. 



I. INTRODUCTION 

Many years ago Weinberg [l| noticed the difficulty to 
account for the 77-77' mass difference. He baptized this 
difficulty as the U(l) problem. Here we address the U(l) 
problem in the microscopic perspective of quark models. 
Notice that all effective models of QCD faced this prob- 
lem. In their second paper, Nambu and Jona-Lasinio 
include the isospin-dependent pseudoscalar contact inter- 
actions in their model to account for the 77-77' mass dif- 
ference, which are reminiscent of scalar and pseudoscalar 
meson exchanges. In this sense it is similar to the origi- 
nal a model of Gell-Mann and Levy [H , where the pseu- 
doscalar exchange does not include the isosinglet. A dif- 
ferent determinant interaction was invented by 't Hooft 
4] to split the 77 and 77' masses. With three flavors this 
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FIG. 1: We show the ladder, minimal geometric series of 
diagrams to include the mesons as Bethe-Salpeter quark- 
antiquark bound states. This ladder is isopin invariant, i.e. 
in the chiral limit the spectrum is identical for the tt, the 
K, the 77 and the 77'. Notice that, starting in this figure, for 
simplicity we depict in the same way the bare and dressed 
gluon-quark-quark vertices. 
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FIG. 2: We show the minimal figure necessary to split the 77 
mass from the 77' mass. This diagram, or any other diagram 
where the in fermion lines are not continued to the out fermion 
lines, only contributes to I = mesons. For other mesons it 
cancels. 



is a tree-body interaction. Recently a four-body interac- 
tion was further added by Osipov, Hiller and Providencia 
H to stabilize the vacuum. At a more fundamental level 
is Lattice QCD, where both the quark and gluon fields 
are considered. However the Lattice QCD simulations 
only reproduce the 77-77' experimental mass difference 
with the indirect technique [6|, |7| of the hairpin diagram 
8[ , or with the extrapolation of the current quark mass 
9(. we aim to understand how the 77-77' mass difference 
arises in a microscopic Quark Model . 

In the Quark Model, and in other related frameworks 
such as the Schwinger-Dyson equations, the minimal 
computations of the meson masses are performed in the 
ladder approximation. The ladder geometric series is de- 
picted in Fig. [TJ In this case the kernel (the set of di- 
agrams that is iterated in the Bethe-Salpeter equation) 
is a one-gluon exchange, or an effective quark-antiquark 
interaction resulting from the integration of the gluon 
field. In the limit of equal quark masses, the ladder ap- 
proximation does not separate the 77 from the 77'. Notice 
that any consistent quark model should also include con- 
sistently chiral symmetry breaking. Then, in the mass 
gap equation, the ladder approximation is equivalent to 
the rainbow, or BCS approximation. In the BCS approx- 
imation, the chiral symmetry is not a SU(Nt) symmetry 
(where Nf is the number of light fermions), but a U(Nf) 



(a) 



(b) 




FIG. 3: In non-perturbative QCD, the minimal diagram 
should be dressed with all possible gluon exchange insertions. 
We show in (a) the minimal one-ladder-figure necessary to 
split the rj mass from the rj' mass. In (b) we show a two-ladder 
diagram, equivalent to the coupling to the initial one-meson 
channel to the two-meson channels. Both these diagrams only 
contribute to / = mesons. 



global symmetry. 

To address the rj - rj' mass difference problem one has 
to go beyond the BCS approximation [HI El]. ^ is neces- 
sary to add at least one quark loop in the mass gap equa- 
tion, including quark-antiquark annihilation and creation 
in the quark-antiquark boundstate equation. The mini- 
mal extension to the kernel of the boundstate equation 
is shown in Fig. [5J Notice that this diagram includes 
two AW trian gles , and thus includes the Adlcr-Bcll- 
Jackiw anomaly |12|. [T3 . [l4| , in this case the non-abelian 
anomaly, also related to the U(l) problem [15]. The 
anomaly is an ultraviolet effect and it is a relevant ef- 
fect. For instance for 77, the anomalous electromagnetic 
decay to 77 is of the same order as the hadronic decays. 

There has been an attempt to compute the diagram 
of FigfT] in SDE approach and explore its effect on 77 - 77' 
mass difference [lg]. However, certain assumptions were 
made about the infrared behavior of the gluon propaga- 
tor, which have been shown later not to be true [171 ] - 

Moreover, in non-perturbative infrared QCD, the dia- 
gram of Fig. [5] is not really complete. The kernel should 
include all diagrams of the same class [l(| • In particular 
the 77-77' mass difference is a non-perturbative problem 
and there is no reason to include only two gluons, or 
two effective quark-antiquark interactions. Inasmuch as 
the full geometric series included in the boundstate min- 
imal boundstate study, the diagram of Fig. [5] should be 
dressed with at least a full ladder series. In Fig. [3J (a) 
we include all the possible number of gluon exchanges, 
and this is equivalent to include a full meson-like ladder. 
This diagram includes all possible t-channel exchanges of 
mesons. 

A different possible way to dress the diagram is de- 
picted in Fig. [3] (b). In this case the diagrams are re- 
summed in two ladders. This is equivalent to couple the 
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FIG. 4: It is also possible to dress the minimal diagram with 
a pure-giuon ladder. We show a different dressing, equivalent 
to couple the pseudoscalar mesons to glueballs. Notice that 
with transverse gluons at least three gluons are needed to 



constitute a pseudoscalar J 



PC 



glueball. This diagram 



also contributes to I = mesons, but it is suppressed by the 
large three-gluon glueball masses. 



original one meson to all possible two-meson channels. 
Notice that, since a pseudoscalar cannot couple to two 
pseudoscalars, the lightest coupled channel includes one 
pseudoscalar and one vector, or one pseudoscalar and one 
scalar. So this channel is already quite heavier than the 
original one meson channel of Fig. [3] (a). In Fig. 2] we 
exchange all the possible number of gluons, not between 
the quark lines, but between the gluon lines, include a 
full glueball- like ladder. This diagram is equivalent to 
couple the pseudoscalar mesons to glueballs. Notice that 
with transverse gluons at least three gluons are needed to 
constitute a pseudoscalar J PC = h glueball. Both in 
lattice QCD and in constituent quark models, including 
models where the gluon mass is generated with a mass 
gap or Schwinger-Dyson equation, the three gluon glue- 
balls arc quite heavy, much heavier than the mesonic cou- 
pled channel of one pseudoscalar and one vector. Thus 
we estimate that the leading non-perturbative diagram 
contributing to the 77 — 77' mass difference is the t-channel 
one meson exchange represented in Fig. [3J (a). In partic- 
ular, because the pion is a light meson, the exchange of 
a virtual pion should be relevant. This is consistent with 
the results of the Nambu and Jona-Lasinio model and of 
the a model. 

Importantly, chiral symmetry forces us to use a self- 
consistent set of diagrams. Similarly to the sets of di- 
agrams that preserve gauge invariance when a photon 
is coupled, in the study of pseudoscalar mesons a self- 
consistent set of diagrams is necessary. 

Notice that different frameworks may be used to com- 
pute the necessary diagrams. We may address the U(l) 
problem in the equal-time quark model formalism, or we 
may use the related euclidian-time Schwinger-Dyson for- 
malism. Although the quark model is explicitly confin- 
ing, and able to reproduce the hadronic spectra up to 
high excitations, the separation of the quark and an- 
tiquark propagators would force us to use a much 
too large number of diagrams. Thus we will work in the 
Schwinger-Dyson formalism, with Euclidian momenta in- 
tegrations, expecting that our results will also apply to 
the quark model. 
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FIG. 5: We show the ladder Ward Identity, crucial to verify 
that the Schwinger Dyson mass gap equation is consistent 
with the Bethe Salpeter bound state equation. 



In Section II we determine the necessary self-consistent 
set of diagrams both for the Schwinger-Dyson quark mass 
gap equation and for the Bethe-Salpeter quark-antiquark 
boundstate equation. Because the full set of diagrams 
goes beyond the present state of the art techniques, in 
Section III we identify the dominant diagrams. Notice 
that the effect of meson exchange in the mass gap equa- 
tion has already been estimated, and thus we focus here 
on the boundstate equation only. We also detail the 
Bethe Salpeter equation, adding the dominant new U(l) 
splitting diagram to the BCS kernel. In Section IV we 
discuss the relevant numerical techniques and we show 
the results for the rj and rj masses. The ideal case of two 
light flavors, relevant to lattice QCD is also explored, 
together with the case of realistic current quark masses 
with three flavors. In section V we conclude. 



II. A CHIRALLY SELF-CONSISTENT CLASS 
OF DIAGRAMS 

We choose the Schwinger-Dyson equation (SDE) for- 
malism to investigate the effect of the diagrams that con- 
tribute differently to the masses of isoscalars (77 or 77') and 
isovector (ir°). For recent reviews on the SDEs and their 
use in hadron physics, see for example [la. Il9j|. 

The SDEs are an infinite number of coupled integral 
equations; a simultaneous, self-consistent solution of the 
complete set is equivalent to a solution of the theory. 
In practice, the complete solution of SDEs is not possi- 
ble for QCD. Therefore one employs a truncation scheme 
by solving only the equations important to the problem 
under consideration and making assumptions for the so- 
lutions of the other equations. 

The simplest Schwinger-Dyson equation is the gap 
equation for the quark propagator. It provides the rela- 
tionship between the quark propagator, the gluon propa- 
gator and the quark-gluon vertex. The exact form of this 
equation is, 



S(p) 1 = Z 2 (ry -p + m bm ) + 



(1) 



j^i9 2 D^{q)—l tJl S{p + q)T a v {q,p), 




FIG. 6: Contribution of the meson exchange kernel to the 
mass gap Schwinger-Dyson equation 

where S{j>) is the flavor dependent fully dressed quark 
propagator, which has the form: 



S(p) 



1 



iM(P 2 ) + B{ P 2 



(2) 



Dfivip — q) is the gluon propagator, and T®(q,p) is the 
quark-gluon vertex, a denotes the flavor of the quark. 
The propagators and the vertex in this equation are 
dressed and renormalized. Z\ and Z 2 are the renormal- 
ization constants of the quark-gluon vertex and the quark 
wave function, and m& m is the bare quark mass. 

The Bethe-Salpeter equation is used to study two- 
particle bound states (mesons). By solving the Bethe- 
Salpeter equation one obtains the Bethe-Salpeter ampli- 
tude which, after proper normalization, completely de- 
scribes the meson as the bound state of a quark and an 
antiquark. The Bethe-Salpeter amplitude corresponds 
to the amputated one-particle-irreducible quark-meson 
vertex, and its Lorentz structure depends on the quan- 
tum numbers of the meson of interest. The exact Bethe- 
Salpeter equation for the meson can be written as, 



d A q 



(2tt)< 



K (q) S(pi + q)T M (pi + q,p 2 + q) S(p 2 



(3) 



where Tjm(j>;P) is the Bethe-Salpeter amplitude for the 
meson M. The Bethe-Salpeter amplitude depends on the 
quark momenta p\ and p 2 . If m is the bound state mass 
then equation ((3|) is only valid for P 2 = (pi— p 2 ) 2 = —m 2 , 
where P is the total momentum of the meson. The kernel 
K{q) is the irreducible quark-antiquark scattering kernel. 

The truncated mass gap SDE and the BSE have to 
be consistent with each other to preserve the chiral sym- 
metry of QCD. This self-consistency is crucial for the 
pion to be a massless boundstate in the chiral limit. It 
has been shown that if one inserts the Bethe Salpeter 
vertex V in all possible dressed quark propagators S(p) 
present in the SDE for the quark self-energy, one must 
recover the kernel of the BSE. Inversely, if one removes 
the Bethe-Salpeter vertex from the BSE, replacing SYS 
by the propagator S, one must recover the self-energy 
of the SDE. This chiral self-consistent relation has al- 
ready been applied to self-energies and to Bethe Salpeter 
kernels containing the infinite ladder series of diagrams 

in hi mi. 

Let us first show that using the rainbow-ladder trun- 
cation (BCS approximation) one obtains a self-consistent 
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FIG. 8: Non-perturbative part of the kernel contributing only 
to flavor singlet rf mass. 



FIG. 7: Contribution of the meson exchange kernel to the 
boundstate Bethe-Salpeter equation 



set of equations for quark propagator, bound states and 
vertices. After that, an analogous set of equations will 
be derived for a different truncation, which includes the 
diagrams necessary for generating 77-77' mass. 

The rainbow truncation of Schwingcr-Dyson equation 
for the quark propagator includes the replacement of the 
dressed quark-gluon vertex by the bare one, 



— — T ^W^-lu, (4) 



where G(q 2 ) is the effective running coupling, q is the 
gluon momentum, and T^ v (q) = — q^q^/q 2 is the 
transverse operator in Landau gauge. The quantity 
T flv {q)/q 2 is the free gluon propagator. 

In the rainbow truncation the equation ([T]) becomes, 

(5) 

where So(p) = (i p + m) is the bare quark propagator, 
m is the current quark mass. 

The equation (|5|) increases the gap between the quark 
and antiquark dispersion relations (which is almost zero 
for bare propagators) and generates the dynamical mass 
M = y 'B 2 (p 2 ) /A 2 (p 2 ) . Because of the color structure 
of the eqn. ((5|), the tadpole does not contribute. Multi- 
plying both sides of eqn. ([5]) with 75 and summing leads 
to, 

S (7Ji) _1 75 +j5S{p 2 y 1 = So(pi) _1 75 +75<S'o(j>2)" 1 



d A q T^ v {q) 
(2tt) 4 q 2 



G {<l 2 )lv [S(pi + 9)75 + 75 S 1 (P2 + q)] lv, 



which is the Bethe Salpeter equation for the vertex, 



Ti(pi,P2) 



{2-Kf 



= 1a(Pi,P2) + 

K (q) Sip, + q)T A (p! 



q,p 2 + q)S(p 2 



(6) 



if the bare 71 and dressed Ta vertices are respectively de- 
fined with the same Axial Ward Identity for an isovector, 

Fa(Pi,P2) = S' _1 (pi)75 +75<S' _1 (P2) , (7) 
1a(jPi,P2) = S , (T 1 (Pi)75 J rlhS^ 1 {p 2 ) , (8) 

and the kernel of the BSE is approximated by the one- 
gluon exchange (ladder truncation, or BCS approxima- 
tion), 

G(q 2 ) A c A c 

K il) -> ^T^(q) Y ^ ® - lv . (9) 

The fact that Axial Ward identity is consistent with 
the ladder approximation for the bound state and the 
rainbow approximation for the quark self energy equation 
ensures that this truncation respects the chiral symme- 
try of QCD. Both approximations are equivalent to the 
planar diagram expansion which is characteristic of the 
Quark Model. 

The bare vertex 7^4 can be computed from the bare 
quark propagator using 

j A (P) = (i / P + 2m)j 5 , P = Pl -p 2 , (10) 

where m is the current quark mass, ja is the particular 
part of the Bethe Salpeter equation for the vertex ([6]), 
and it vanishes when the current quark mass m is small 
(chiral limit) and at the same time the total momentum 
P M of the vertex is small. 

At this point it is important to clarify that in general, 
as it follows from (fTU|) . the bare vertex ja is a combina- 
tion of the axial vector vertex 7^75 and the pseudoscalar 
vertex 75. In the chiral limit 7^4 becomes pure axial vec- 
tor vertex contracted with momentum, and in the limit 
of vanishing momentum P M it becomes pure pseudoscalar 
vertex multiplied by the current quark mass. For simplic- 
ity (and because they are defined with the Axial Ward 
identity ([7]) and (HJ)), in the rest of this paper 7^4 will 
be called the bare axial vertex and Fa will be called the 
dressed axial vertex, although they possess a more gen- 
eral Dirac structure. 
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FIG. 9: Non-perturbative part of the kernel shifting all the 
pseudo-scalar masses equally. 

The dressed vertex Ta can be computed from the 
dressed quark propagator @ using 0, 

T A (p u p 2 ) = [iA{p 1 )i> 1 -iA{p 2 )i> 2 + B{p 1 ) + B{p 2 )] lb . 

Ta is finite providing spontaneous chiral symmetry 
breaking occurs in eqn. ([5]) and a dynamical mass of 
the quark is generated. For instance, if the total momen- 
tum P = pi — p 2 of the vertex vanishes, the vertex Ta is 
simply identical to 25(75)75, where B{p) is a finite solu- 
tion of the mass gap equation. This vertex describes the 
coupling of the pion to two quarks, and eqn. ([6|) shows 
that in the chiral limit we have a Goldstone boson, the 
pion. 

For simplicity the flavor is not yet included. Flavor will 
only be explicitly included at the end of subsection VII. 
The isoscalar axial Ward identity must include the Axial 
anomaly, which is crucial to the U(l) problem. Never- 
theless the pion is an isovector, and in the coupling of a 
pion we do not need to concern with the Axial anomaly. 

We now derive a powerful Ward identity for the ladder 
which involve the axial vertices and the ladder. This 
relation is depicted in Fig. [5l it constitutes an extension 
of the Ward Identity for the propagators in eqs. Q and 
([5]). This identity is derived if we expand [hJ 0, HH 
the ladders and substitute the vertex in the left hand 
side. Then all terms with an intermediate 75 include 
the anticommutator {75 , 7 V } and this cancels because 
the interaction is chiral invariant and the kernel is local. 
Only the right hand side survives. 

We are now ready to derive the mass gap and the 
boundstate equations beyond rainbow-ladder truncation, 
to include the diagrams that are able to generate 77-77' 
mass difference. Starting from the diagram in Fig. [3] 
(a) (which is the minimal diagram that involves one t- 
channel meson exchange), and substituting the Bethc- 
Salpeter vertex by a propagator, we arrive at the contri- 
bution for the quark self-energy in Fig. [5] Notice that, 




FIG. 10: The homogeneous Bethe-Salpeter equation we solve 
here. The kernel includes the BCS level One Gluon Exchange, 
the anomalous annihilation into two gluons, and the non- 
perturbative One Pion Exchange. 



to remove any disconnected diagrams, the first diagram 
must include at least one gluon exchange, thus, in Fig. 
[6l the geometric series starts with one gluon exchange. 

Then, inserting the Bethe-Salpeter vertex in all pos- 
sible propagators of the quark self-energy, we get the 
kernel diagrams in Fig. [7] for the boundstate equations. 
The first kernel diagram is obtained inserting the vertex 
in the only propagator of the self-energy exterior to the 
ladder. The second kernel diagram is obtained inserting 
the vertex in the same quark line, but inside the ladder 
of the self-energy. The third kernel diagram is obtained 
inserting the vertex in the quark line linking to the ex- 
ternal legs of the self-energy. We can also express these 
diagrams in terms of ladders only, and we finally get the 
diagrams in Figs. [5] and [5] In Fig. [5] we show the di- 
agrams contributing to the U(1)a mass splitting and in 
Fig. [9] we show the diagrams contributing to all mesons. 



III. A CALCULABLE BETHE SALPETER 
KERNEL 

We now choose the best framework to compute the new 
beyond BCS diagrams. Notice that the diagrams in Figs. 
[6l [8] and [9] all include full ladders, and an integral in one 
of the relative variables of the full ladder. Technically 
this remains a problem since it goes beyond the present 
state of the art of quark models. 

Two different classes of quark models, consistent with 
chiral symmetry, have been investigated and applied to 
different problems. The equal time quark models are 
confining, thus they are adequate for the computation 
of the full ladder, including the full spectrum of mesons. 
However they are not Lorentz invariant and thus are in- 
adequate for the boost of the ladder, occurring in the mo- 
mentum space integral in the Feynman diagrams. And 
the Lorentz non-invariance also splits the time-like from 
the space-like components and observables of the model. 
On the other hand, the euclidian space quark models 
are convenient for full integrations of internal momenta 
in Feynman diagrams, but are not confining and thus 
are inadequate for the computation of the full ladder. 
However, the full ladders can be replaced by the lowest 
boundstate contribution, i. e. approximately assumed to 




FIG. 11: An example of BSE equation needed to solve for 
the mass of the bound state uu. In the annihilation dia- 
gram quark propagators of all flavors contribute. That makes 
it necessary to include kaon exchange in this diagram if the 
strange flavor is taken into account. 
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be similar to the pseudoscalar pole contribution. In this 
case the euclidian models, essentially adequate to address 
the lowest energy phenomena of hadronic physics (includ- 
ing the pseudoscalar ground states to the vector ground 
states and pseudoscalar first excitations), can be applied 
to the U(1)a breaking of the pseudoscalar spectrum. 

The dominant part of the non-perturbative dressing 
of the kernel of the Bethe Salpeter equation with hadron 
ladders is the One Meson Exchange in the t-channel. No- 
tice that the corresponding diagram occurs with two mi- 
nus signs when compared with the One Gluon Exchange. 
A first minus sign is necessary to cancel all the discon- 
nected diagrams. A second minus sign is due to the 
fermion loop. 

Now, t-channel exchange is dominated by the pseu- 
doscalar meson exchange, since the pseudoscalar mesons 
are the lightest ones. Thus, in this first attempt of the rj 
mass study beyond BCS, we will approximately saturate 
the meson exchange with only the ground state pseu- 
doscalar nonet. And, because we are mostly interested 
in the U(1)a breaking, we will also neglect all terms that 
contribute equally to all mesons, providing only a con- 
stant overall mass shift of the different spectra. Thus 
we will not consider the contribution to the mass gap 
equation of Fig. [5] and the contribution to the bound- 
state kernel in Fig [9] (it has also been estimated by one 
of us [l0| that these effects are relatively small) . In any 
case, these two effects cancel in the chiral limit, since 
our self-consistent approach complies with the Goldstonc 
theorem. 

The minimal kernel necessary to estimate the differ- 
ent contributions to the 77 and 77' masses is presented in 
Fig. [TU] The first diagram on the right-hand-side of Fig. 
[10] is the BCS level One Gluon Exchange. This diagram 
contributes to all mesons. The second diagram on the 
right-hand-side includes the pion exchange in t-channel 
between quark and antiquark. It only contributes to the 
isosinglet mesons, and therefore generates the mass dif- 
ference between the ttq, the 77 and the 77'. 

The diagrammatic form of Bethe-Salpeter equation de- 
picted in Fig. QJj]can be written as, 



FIG. 12: Singularities of the solution of the gap equation 
for u quark. The real part of the scalar function a v (p 2 ) is 
presented, which is defined in (I16|l . 



(2tt)< 



K (p, q- P) S a (q + V P)T a v a (q; P) S a (q - 0P) 



where 

Gt (p, 9; P) = r#(p, q, P)G M (p - q)T%(p, q, P), (13) 

r^J is the Bethe-Salpeter amplitude and Gm is the prop- 
agator of the meson being exchanged, a and b denote 
quark flavors. 

The effective running coupling in the kernel can be 
modeled to account for the effects of the truncation. We 
used the Maris-Tandy model [22, 23] in our calculations, 



G(k 2 
k 2 



4n 2 D 



4ir 2 lm F(k 2 ) 



I In 



k 2 /A 2 



QCD 



(14) 
(15) 



1 v 



(12) 



where F{k 2 ) = (l - exp^j /k 2 , 7m = 12/ (33 - 2N f ) 

an where the parameters are m t — 0.5 GeV, Nf = 4, 
Aqcd = 0.234 GeV, r = e 2 - 1, D = 0.93 GeV 2 and 
w = 0.4 GeV. This model has been shown to work well 
for the description of the light meson properties (see for 
example 0, HH|). Therefore we expect it to give sensible 
results for the masses of 77 and 77'. 

In general, for Nf quark flavors there will be Nf cou- 
pled integral Bethe-Salpeter equations, which have to be 
solved together. An example of such an equation for the 
uu bound state with all possible pseudoscalar meson ex- 
changes is depicted in Fig. [TTJ 

Notice that one of the features of Schwinger-Dyson 
equation approach is the fact that the solution for the 
quark propagator has singularities. An example of the 
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FIG. 13: Location of the singularities of the quark propaga- 
tor with respect to our largest integration domain of Bethe- 
Salpeter equations. 



solution for the function a v (p 2 ) for the u quark, 



A(p 2 ) 



p 2 A 2 {p 2 ) + B 2 {p 2 )' 



(16) 



is presented in Fig. 1121 In general, these singularities can 
pose a serious technical problem for the numerical calcu- 
lations if present in the integration domain of the inte- 
gral Bethe-Salpeter equations. However, we have checked 
that in our case these singularities are outside from the 
largest integration domain that we have. Our largest in- 
tegration domain and the positions of the poles are shown 
in Fig. [13] 

Thus the masses of the rj and rj mesons are computable 
in the SDE and BSE formalism, extending the Klabucar 
and Kekez parametrization of the U a (l) breaking 27]. 



IV. RESULTS 

We solve the Bethe-Salpeter equation ([12")) for the pseu- 
doscalar mesons. This equation takes into account the 
additional "annihilation" diagram which contributes to 
the isoscalar but not to the isovector mesons. 

This additional diagram is depicted in Fig. [3^i, it con- 
sists of the annihilation of our qq pair into an infinite 
gluon ladder. If calculated perturbatively, this diagram 
does not contribute to the isoscalar meson mass in the 
chiral limit. However, it gives rise to the rj' meson mass 
if treated nonperturbatively. We propose to represent 
the infinite gluon ladder in the "annihilation" diagram 
as a t-channel meson exchange. In general, mesons in a 
t-channel can have any quantum numbers. However, the 
pion and the other members of the pseudoscalar octet, 
being the lightest mesons, give the dominant contribu- 
tion. 




FIG. 14: Dependence of the mass of rj meson on the "weight" 
of the additional diagram in the chiral limit for two and three 
quark flavors. 



This Ua{1) breaking term is added to the ladder Bethe- 
Salpeter kernel in the Bethe-Salpeter equation (TT2"]) , as in 
Fig. [TUJ The Bethe-Salpeter amplitude, needed for the 
pseudoscalar-exchange kernel, is previously calculated by 
solving Bethe-Salpeter equation in the ladder truncation 
©■ The quark propagator, which is needed for the solu- 
tion of the equation (fT2|) is also obtained by solving the 
ladder Schwinger-Dyson equation {1} numerically. This 
ensures that the pion and kaon non-isoscalar masses are 
unaffected. Having thus defined both the Bethe-Salpeter 
kernel and the quark propagator, we solve numerically 
the Bethe-Salpeter equation by the power method. 

First, we investigate the isoscalar meson mass in the 
chiral limit. The chiral limit is the limit where the cur- 
rent quark masses vanish. In this case the " annihilation" 
diagram is simply proportional to the number of flavors 
Nf. From the extent of the J7a(1) breaking it is expected, 
that when the chiral limit is taken for all three quark fla- 
vors, the 7To and rj should become massless, and the rj' 
should remain massive. We are able to reproduce this be- 
havior with our phenomenological model: the " annihila- 
tion" diagram does not contribute to the ttq and rj in the 
chiral limit leaving these mesons massless, but changes 
the mass of the flavor scalar rj' (Fig. . 

The Bethe Salpeter amplitude of pseudoscalar 
mesons can be decomposed into four Dirac structures 
E, F, G, H. By keeping only the E amplitude (separa- 
ble in the chiral limit and dominant for realistic masses) 
we calculate the rj' meson mass in the chiral limit with 
the dependence on the "weight" of the additional dia- 
gram (phenomenological factor / which multiplies the 
additional U a{\) breaking diagram in the BSE). The re- 
sult for 2 and 3 quark flavors is presented in Fig. 03] In 
the two-flavor case there is only one rj, in Fig. Q3] it is 
also denoted rj '. One can see that our calculations give 
a very reasonable estimate of the mass of rj in the chiral 
limit. The solutions of the Bethe-Salpeter equation for 
7T, K and rj' (in the chiral limit) are presented in Fig. 
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FIG. 15: The E(p 2 ) Bethe-Salpeter amplitudes, solutions of 
Bethe-Salpeter equation, for the n, the K and the rj' . 
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FIG. 16: Dependence of the flavor scalar r/o mass on the 
"weight" factor / for two and three quark flavors. 



[T31 Notice that in the chiral limit the pion and kaon 
amplitudes are identical. 

Next, in our calculation, we employ broken SU(3)f 
symmetry with u and d quarks of equal finite mass, and 
realistically heavier s quark. That means that in the " an- 
nihilation" diagram we have to include both kaon and 
pion propagators. We examine the effect that this dia- 
gram has on the physical mesons. It has been shown be- 
fore [24[ that the SDE and BSE approach, in the rainbow- 
ladder truncation, works reasonably well for the prop- 
erties of the light isovector mesons, and we reproduce 
these results obtaining for the pion = 0.139 GeV and 
_U = 0.131 GeV. 

The SU(3)f octet and singlet isospin zero states, 77s 
and 770, can be expressed in the qq-basis, 

\m) = 4= {\uu) + \dd) - 2\ss)) , (17) 
\m) = (M + \dd) + \ssj) . (18) 

while in SU(2)f only the 770 state is defined. The "anni- 
hilation diagram" does not contribute to the mass of 77s 
in the chiral limit. For the finite quark masses, however, 
it does still contribute due to the mass difference of u(d) 
and s quarks. For 770 this diagram makes a difference 
both in chiral limit and in the finite quark mass case. To 
find the masses of the rjo and rjg we have to solve the 
system of two coupled integral BSEs (one for u(d) and 
one for s flavor), equivalent to the equation (fT2f . An ex- 
ample of such an equation for the uu channel, including 
its coupling to the dd and ss channels, is depicted dia- 
grammatically in Fig. 1111 The dependencies of 770 and 
77s meson masses on the "weight" factor / are presented 
in Fig. [TBI and [T71 The full BSE has been solved for each 
"weight" factor. 

To make a connection with the physical mass eigen- 
states, we introduce the mixing angle O, 



\V 



sin#|?7 8 ) + cos9\i]o) 



(20) 



By comparing the results of our calculation to the ex- 
perimental values of 77 and 77' mass we can determine the 
mixing angle 



tan 2 9= M ^~ M * 
M v , - M m 

tan 2 e- M "'- M - 



M, 



■t]u 



(21) 
(22) 



\v) 



cos6/|7y 8 ) - sin 1 770) 



(19) 



For the "weight" factor / = 0.9 we obtain 9 » -28° in 
reasonable agreement with experiment, which favors the 
mixing angle in the vicinity of —20° (27l.[28l|. Now we can 
use this angle to explore the dependence of the physical 
rj and 77' on the weight factor. The corresponding results 
are presented in FigfT8l 



V. CONCLUSION 

We identify a full one-quark-loop class of microscopic 
quark diagrams contributing to the 77 and rj' masses, sep- 
arating these pseudoscalar mesons from the ir and K 
mesons. 

Because it is not yet possible to compute the full class 
of diagrams, we identify the dominant diagrams, equiv- 
alent to the exchange of the ground state pseudoscalar 
mesons. 

Nevertheless this problem still requires state of the 
art techniques. We avoid the equal-time framework, 
where the number of diagrams would be much larger be- 
cause the quark propagators would be separated from 
the antiquark propagators. We adopt the Euclidian time 
Schwinger Dyson framework. The Bethe-Salpeter equa- 
tion turns out to be solvable, because the poles in the 
solution of the gap equation, for negative Euclidian mo- 
mentum, does not cross the complex domain of integra- 
tion in the Argand plot of the momentum, depicted in 
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FIG. 17: Dependence of the r/s mass on the "weight" factor / 
for three quark flavors. 

Fig. [13] Together with the isovector ground state pseu- 
doscalars [23|,[24],[25[, first excited pseudoscalar states [IB] 
and ground state vectors (22j, this is a new case among 
the few where the Euclidian Bethe Salpeter equation has 
been solved so far. 

We notice that the infrared part of the complete class 
of diagrams should vanish in the chiral limit. Only the 
ultraviolet anomaly of Adler Bell and Jackiw should sur- 
vive in that limit. Thus it is a priori expected that our 
dominant diagram should be excessive to produce the 
r\ and rj masses, since a cancelation with the other dia- 
grams should occur in the chiral limit. Hence we multiply 
this diagram by a reducing factor /. We find that our 
diagram is indeed dominant since only the factor / = 0.9 
is necessary to arrive at the correct experimental masses. 
Importantly, we then essentially comply with the exper- 
imental mixing of the 770 and ry 8 . 

Our results show that the quark model may solve the 
U(l) problem with simple microscopic interactions, al- 
though new technical advances in the computation of full 
ladders are necessary before this problem is fully solved. 
More progresses may be achieved with the computation 
of contributions to the diagrams in Fig. [3] (b) and with 
the separation of the infrared and of the ultraviolet con- 
tributions to the 77 and rj masses. 



Importantly we also study the case where there are 
only two light flavors. This case is not realistic, never- 
theless it is important to be compared with the two-flavor 
studied in Lattice QCD. In this two-flavor case we find 
that the 77' mass is only reduced by 20% . This stresses 
the relevance of two-flavor studies in Lattice QCD. 

On the other hand, the importance of the one-pion 
exchange diagram explains why Lattice QCD, with the 
present size of the lattices, cannot directly reproduce the 
77 and 7/ masses. Since the pion is very light, a large 
lattice will be necessary to encompass with no chiral ex- 
trapolation the extent of the one pion exchange potential 
between the quarks. 
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